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A TECHNIQUE  FOR  SONAR  DATA  PROCESSING 


by 


W.  J.  Dejka 


PREFACE 


This  memorandum  describes  a novel  approach  to  sonar 
data  processing  growing  out  of  study  by  the  author  of 
the  sonar  requirements  assumed  for  the  ASM  Hydrofoil 
during  work  on  the  command  and  control  subsystems  for 
this  craft  type.  This  memorandum  has  been  prepared  be- 


cause It  Is  felt  that  the  Information  may  be  of  use  to 
other  Laboratory  employees  In  Its  present  form,  ^yery 


- pie  LuffT 


This  memorandum  should  not  be  construed  as  a 


report,  as  Its  only  function  Is  to  present  for  the  Infor- 


mation of  others  a small  portion  of  the  work  being  done 
on  this  problem  (HEL  Problem  j6-l;  BU5HIPS  SS  600  000, 
Task  1723). 


INTRODUCTION 


\ 

In  the  command  and  control  portion  of  the  ASW  Hydrofoil 
feasibility  study  fMEL  Problem /j6«-l^^equirement8  for  sonar 
detection  and  tracking  were  postvilated  (-NEL^EesearGh  Reports  '' 
-1118,— The  problem  of  extracting  information  from  the 
sonar  signal  returns  for  tracking  purposes  were  intensified 
by  the  high  target  closing  rates,  the  modest  sonar  ranges  and 
attendant  high  ping  rates,  and  the  need  to  provide  multiple- 
target-tracking  capability.  The  data  processing  rate  took 
incremental  steps  as  these  requirements  took  form.  Although 
certain  aspects  of  the  LORAD,  SPADE  and  ASIAC  developmental 
sonar  data  processing  systems  are  probably  applicable,  as 
further  study  may  show,  none  of  these  as  a system  fully  sat- 
isfy the  ASM  Hydrofoil  sonar  data  processing  requirements. 

It  has  been  established  that  a coherent  detector 
(matched  filter)  will  be  used  in  the  sonar  system,  and  its 
characteristics  are  well  known.  Consequently,  beginning  with 
its  output  the  subsequent  sonar  data  processing  specifica- 
tions are  to  be  determined,  as  well  as  the  design  approach 
required  to  meet  these  specifications,  inie  following  discus- 
sion will  atten^t  to  give  some  insight  into  this  particular 
problem. 


Problem 


The  design  of  a sonar  data  processor  specifically  to 
perform  the  functions  of  detection,  classification  and 
tracking  is  of  present  concern  in  the  ASW  Hydrofoil  feasi- 
bility study.  Though  present  developmental  processing 
systems  have  been  somewhat  successful,  none  is  directly  ap- 
plicable as  a solution  to  the  AHC  sonar  processing  problem. 
Some  specific  as8Uii5)tions  which  prcmqpt  new  approaches  to 
processing  development  are: 

1.  The  data  (ping)  rate  of  the  sonar  is  extremely 
high  (range  is  relatively  short,  about  3000  yards). 

2.  Detection  must  be  automatic  or  semi-automatic  and 
it  is  to  be  completed  in  three  pings. 

3.  Classification  is  to  be  completed  in  five  pings  (or 
in  20  secs  at  a 14-  sec.  ping  rate). 

4.  Tracking  accuracy  for  fire  control  use  will  have 
reached  its  maximum  in  twelve  pings  (48  secs  until  a weapon 
can  be  fired). 

5.  The  advanced  hydrofoil  craft  will  be  space  limited, 
and  the  techniques  mentioned  above  suffer  from  ccxaplexity 
and  size  limitations. 


2 


The  basic  approach  that  will  be  taken  here  is  as  follows:  . 


A model  of  the  target,  environment,  and  soneir  are  to  be  devel- 
oped which  will  relate  all  observations  (or  clues)  into  one 
model  of  a randcm  process.  This  model  will  be  a time  depend- 
ent reuidom  process,  that  is,  seme,  or  possibly  all,  observa- 
tions will  be  a function  of  time.  The  model  will  be  non-linear 
but  can  be  linearized  for  purposes  of  simplifying  computation. 

This  model  will  then  be  employed  in  ceLlciilating  the  maxi- 
mum-likelihood ratio  for  best  estimates  of  the  sonar  data  or 
information  that  can  be  extracted  from  the  clues  or  observations. 
It  will  be  developed  in  a manner  which  gives  a recursive 
conqmter  algorithm. 

Statistical  Estimation  of  Parameters — Likelihood  Function 


The  theory  of  the  estimation  of  parameters  is  part  of  the 
mathematical  theory  of  statistics  and  is  utilized  frequently 
in  the  communications  and  information  processing  fields.  In 
sonar  data  processing,  a number  of  observations are 
available  at  discrete  time  increments.  These  time  dependent 
observations  are  clues  which  must  be  interpreted  to  evaluate 
whether  a target  is  sub  or  non-sub,  and  whether  it  is  hostile 
or  friendly.  This  interpretation  of  parameters  frexn  given 
observations  can  be  treated  by  the  method  of  statistical 


estimation  of  time-dependent  parameters. 


lIMs  treatment  requires  one  major  a8SUBq>tion  tihlch  x>ennlt8  the 
problem  to  be  solved  by  analog  as  veil  as  digital  equipment.  It  is 
assumed  that  a complete  statistical  description  of  the  error  or  noise 
Is  known.  For  most  purposes^  the  noise  or  error  may  be  assumed  Gaussian. 
Gaussian  noise  Is  representative  of  memy  practical  noise  processes  and^ 

In  addition^  makes  the  mathematical  an^Jysls  less  coopUcated  than  do 
other  distributions.  This  assumption  Is  added  to  asstxre  that^  should 
the  real  time  solution  on  a digital  computer  be  too  time  consuming,  an 
analog  liiq>lementatlon  is  feasible. 


The  observations  are  represented  in  vector-matrix  notations  by 

y(t)  = £1,  

.here  le  hhe  sum  of  the  ooiee  andasignal^ 

The  signal  or  parameters  to  be  estimated  are  related  and  this  functional 
relation  Is  assvoned  linear.  The  observations  are  made  at  discrete 
IntervaDspr  else  continuously  for  a finite  time  Interval.  The  discrete 


form  is  represented  by  ^ ^ ^ TX  £ 


<;■=  4 . . ,-t' 


These  N valxies  of  used  to  make  estimates  of  the  parameters. 

Since  the  complete  statlstlceG.  description  of  the  noise  is  assumed 
known,  the  Joint  probability-density  function  of  noise  . /^h  ^ ^ 

is  available  to  the  observer.  And,  since  ^ t Joint 

probability-density  function  for  the  N observed  veO-ues  of  the 
sighhi  Lj(t)  i.  - 

~ i/n  ^ 

Therefore,  the  function  of  the  estimator  is  to  form  estimates  of  the 
(Ji\  on  the  basis  of  the  N observations  / ** ' described 


k 


whose  form  Is  Imovni  and 


by  a probability-density  function  fs 
depends  on  the  CLi  • 

Let  be  any  estimate  of  a single  parameter  d-j  based  on 


the  N observations 
the  inequality 


T5ie  variance  of  CLj,  obeys 

X id 


Variance  ^ “ 


\_zecL^ 


where  E denotes  the  expected  or  mean  value  of  a random  variable.  If 
the  expected  value  is  equal  to  the  parameter  of  interest  [that  is,E(a)=a 
for  n.n  a's]  the  estimates  are  called  unbiased  and 


Model  of  a Random  Process 

The  development  of  a model  of  a random  process  to  represent  the 
submarine  target  in  a changing  ensonified  environment  heis  been  the 
subject  of  extensive  sonar  studies.  Ibough  an  accurate  model  is  highly 
desirable  for  increased  accuracy,  the  approach  to  follow  allows  a crude 
linear  model  to  be  postulated, and  then  interprets  it  as  a nonstationary 
model,  niis  interpretation  permits  the  model  uncertainties  to  be  attri- 
buted to  this  stochastic,  non-stationary, property. 

Consider  the  following  simple  example.  Tlie  model  is  a first-order 
difference  equation. 
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where  x Is  a parameter  or  variable  of  interest.  The  random  variable, 
is  the  course  of  the  non- stationary  property  of  the  model.  The 
variance  and  mean  of  sure  known  or  can  be  approximated  from 


experimental  work. 

OMs  linear  model  is  quite  versatile  and  ceui  be  used  satisfactorily 
in  most  problems.  Should  experimental  use  of  this  proposed  system  prove 
inadequate,  a non-linear  or  time-varying  process  model  can  be  utilized 
without  undue  complexity. 

Observations  are  given  hy  ^ ~ ^ I? 


such  that  they  are  linear  combination  of  the  process  state  ^ 
measurement  or  equipment  error,  • Here  agedn  the  statistics 

of  these  random  variables  are  experimentally  determined. 

Uiis  model  is  now  used  in  the  statistical  estimator  for  determining 
the  maximum  likelihood  of  the  states  (clues)  for  classification. 

Discrete  Versus  Continuous 

The  computation  time  as  well  as  the  conqiutational  storeige  requirements 
indicate  that  the  conqzlex  system  should  be  both  discrete  and  continuous 
with  a binary  decision  as  an  output.  Such  an  output  is  sub(l)  or  non-sub  (o) 
and  it  should  automatically  activate  an  alarm.  The  fact  that  a hybrid 
system  is  required  is  not  a limitation  but  may  actually  be  an  asset. 

Such  decision  techniques  have  been  sadly  neglected  and  they  hold  much 
promise,  at  least  mathematically. 

A disadvantage  to  modeling  the  random  process  as  a hybrid  system 
is  that  the  system  designers  frequently  are  skilled  in  either  discrete 


k: 


and 
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or  continuous  design  but  not  both.  This  should  not  be  a limitation^ 


for  those  >rtio  attempt  such  hybrid  modeling  are  revrarded  with 


some  Interesting  results.  Practically,  it  may  be  the  only  way  the 


problan  can  be  iiiQ>lemented. 


A Linear  ^stem  as  a Likelihood  Decision  Device 


Having  applied  some  smoothing  to  a set  of  clues  (the  maximum -likeli- 
hood estimates), It  is  necessary  to  evaluate  whether  the  target  is  sub 


or  non-sub  given  the  above  clues  or  observations.  It  is  desired  to 


perfom  this  decision-making  continuously , as  along  a given  sonar  beam. 


To  show  how  this  is  accomplished, consider  the  simplest  hypothesis 
testing  problem.  Let  the  observed  signal  xCt)  be  due  to  noise 


only  or  to  a precisely  known  signal  and  noise.  The  former 


hypothesis,  noise  only,  is  denoted  by  and  the  latter  hypothesis. 


signal  and  noise,  is  denoted  by  H|  ,It  is  desired  to  devise  a test 
for  deciding  in  favor  of  or  H|  , 

If  one  lets  be  the  probability  density  that  if  is 


true,  the  observed  wavefono,  x(t)  could  have  arisen;  and  if  f,(x) 


is  the  probability  density  that  if  )4,  is  true,  x(-t]  could  have 


arisen;  then  the  test  has  the  form  \-)  (x\ 

aceefT  W,  i-f  > A 


Here  X is  a constant^  dependent  on  "a  priori"  probabilities  and  costs. 


if  these  are  known,or  on  the  predetermined  value  of  oC  , the  false 


alarm  probability.  !Bie  test  asks  us  to  examine  the  possible  causes  of 


what  we  have  observed. 


L 
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Let  us  assume  now  that  the  noise  n('^)  is  additive,  gaussian 
and  white  with  spectral  density  ( ) and  further  that  the  signal, 

if  present,  has  the  known  form  ^ where 

the  delay  and  the  signal  diiration  T are  assumed  known.  Ilien,  on 
observing  x.(t)  in  some  observaticn  interval,  X”,  which  includes  the 
Interval  y/T’  / the  two  hypotheses  concerning  its  origin  are: 

Ho  ' ^(^)=  nCir)  ± ifj  I 

hf , : Y-  fr-)  - ^ ) -f  /L  ftJ  I 

Now,  the  probability  density  of  a sample,  , of  white  gaussian 


noise  lasting  from  a to  b may  be  expressed  as 

ij  " 

where  is  the  double-ended  spectral  density  of  the  noise,  and  k is 

a constant  not  dependent  on  ■hM  , Hence  the  likelihood  that,  if 
f/c  ts  true,  the  observation  could  arise  is  simply  the 

probability  density  that  the  noise  waveform  ceinessume  the  form  of 

i.e.  -p^^x)  = 1?  exf  J 


the  region  of  integration  being,  as  indicated,  the  observation  interval,  I. 
Similarly,  the  likelihood  that,  if  /4,  is  true,  could  arise  is 

the  probability  density  that  the  noise  can  assume  the  form 

where  we  have  denoted ^ energy  of  the  signal,  by  E. 

On  substituting  the  latter  expression  and  taking  the  logarithm  of 


both  sides  of  the  inequality,  the  hypothesis-testing  criterion  becomes 


inequality,  the  hj^othesis-testing  ( 

H,  It  cj  ("tc ) > 

accept-  f+o  ;f  cj  (O  X'' 


